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Recent discoveries in neutron scattering experiments for Kapellasite and Herbertsmithite as well 
as theoretical calculations of possible spin liquid phases have revived interest in magnetic phenom- 
ena on the kagome lattice. We study the quantum phase diagram of the S = 1/2 Heisenberg kagome 
model as a function of nearest neighbor coupling J\ and second neighbor coupling J2. Employing 
the pseudofermion functional renormalization group, we find four types of magnetic quantum order 
(q = order, cuboc order, ferromagnetic order, and \f?> x a/3 order) as well as extended magneti- 
cally disordered regions by which we specify the possible parameter regime for Kapellasite. In the 
disordered regime ^ < 1, the flatness of the magnetic susceptibility at the zone boundary which 
is observed for Herbertsmithite can be reconciled with the presence of small J 2 > coupling. In 
particular, we analyze the dimer susceptibilities related to different valence bond crystal (VBC) 
patterns, which are strongly inhomogeneous indicating the rejection of VBC order in the RG flow. 

PACS numbers: 75.10.Jm, 75.10.Kt 



Introduction. Frustrated magnetism is a focus of con- 
temporary research in condensed matter physics, com- 
bining a plethora of experimental scenarios and diverse 
theoretical approaches to describe them. One of the most 
fascinating challenges of the field has been to investigate 
and understand the interplay of magnetic quantum order 
and disorder on the kagome lattice. A major reason why 
this lattice of corner-sharing triangles yields such a com- 
plicated structure of quantum phases is already evident 
from the classical kagome Heisenberg model (KHM): As 
a function of nearest neighbor and next nearest neigh- 
bor Heisenberg couplings J\ and J2, many different mag- 
netic orders are present pQ, where an infinite number of 
degenerate ground states can be found [2] . From a theo- 
retical perspective, not many rigorous results about the 
quantum phase diagram are known so far. Advanced 
mean field theories have provided important guidance as 
to what type of ordered and disordered quantum phases 
could possibly be found [3-8 , but cannot give unam- 
biguous information about which phase will eventually 
be stabilized in the microscopic model. A peculiar fea- 
ture of the KHM which is known since early exact nu- 
merical calculations of finite size clusters [9] is the large 
amount of singlet states at low energy. This suggests 
a plethora of competing quantum-disordered phases and 
is probably one of the main reasons why the interpreta- 
tion of present results for the J\ KHM from microscopic 
numerical approaches is not yet settled [T0HT3] . 

Transferring our fragile theoretical knowledge to exper- 
imental scenarios is even more challenging. The Herbert- 
smithite compound ZnCu3(OH) 6 Cl2 is one of the rare 
properly investigated material realizations of a S = 1/2 
kagome spin model, which is supposedly dominated by 
the Ji term. Early investigations from neutrons [14] and 
muon spin rotation [15] have already indicated its uncon- 



ventional frustrated magnetic properties, exhibiting no 
sign of magnetic order down to a few mK. While there 
is no indication for a magnetic spin gap, it is likely that 
this does not hint at a key property of a possible KHM 
description, but might be due to Dzyaloshinskii-Moriya 
(DM) effects and impurities which further complicate the 
picture [TJJ H6HT9] . Detailed latest neutron scattering 
experiments resolve an extremely fiat magnetic suscep- 
tibility profile [20] [21], which allows for an interpreta- 
tion along spin fractionalization as known from spinon 
continua in the spin structure factor of quasi- Id sys- 
tems [22 . Recently, neutron and /iSR experiments on 
Kapellasite Cu3Zn(OH) 6 Cl2, of which Herbertsmithite 
is a polymorph, have found indication for cuboc mag- 
netic order [23] . This is an exciting discovery, as elec- 
tronic structure calculations find dominant antiferromag- 
netic Ji coupling and hence a similar regime as Herbert- 
smithite [24] . From the knowledge where classical cuboc 
order emerges, however, this suggests the presence of fer- 
romagnetic J\ and antiferromagnetic J2 in Kapellasite. 

In this Letter, we develop a pseudofermion functional 
renormalization group (PFFRG) perspective on the J\- 
J2 kagome Heisenberg model which is ideally suited to 
treat magnetic order and disorder tendencies on an unbi- 
ased footing [25 -29 . Our objectives are two-fold. First, 
we obtain a detailed understanding of the J1-J2 quan- 
tum phase diagram, including all magnetically ordered 
and disordered regimes as well as the associated phase 
transitions. In particular, in light of new findings for 
compounds such as Kapellasite, we allow both J\ and 
J 2 couplings to be ferromagnetic and antiferromagnetic. 
Aside from ferromagnetic, cuboc, y/3 x and q = 
order, we find magnetically disordered regimes located 
around (Ji, J2) ~ (1,0), (Ji, J2) ~ (0, 1), and possibly a 
nonmagnetic phase separating the cuboc from the ferro- 




FIG. 1. Top center: Phase diagram of the J1-J2 kagome model. Colored regions correspond to the magnetically ordered phases 
q = 0, cuboc, ferromagnetic and x V3 order, white regimes are non-magnetic. Faint areas near the phase transitions (e.g. 
between the cuboc and the ferromagnetic phase) are regions of enhanced uncertainty. Bottom center: Flow of the susceptibilities 
in the four ordered regimes. Arrows indicate the instability features. The curves for q = 0, cuboc, ferromagnetic and V3 x y/3 
order are given for a — 45°, 130°, 180° and 315°, respectively (black dots in the phase diagram). The left and right sides of 
the figure show the susceptibility profiles for these types of orders (at the same values for a) at the instability breakdown. 



magnetic domain (Fig. [I]). Second, we specifically investi- 
gate the Ji-KHM disordered regime which is supposed to 
relate to the Herbertsmithite scenario. Due to the large 
system sizes of up to 317 sites which are reached by PF- 
FRG, we obtain accurate resolution of the momentum- 
resolved static magnetic susceptibility. Aside from short- 
range correlations, we observe a broad spectral distribu- 
tion which becomes flat at the magnetic zone boundary 
for J2/J1 - 0.017 (Fig. [2}. This profile is similar to 
what is observed in recent neutron measurements [21]. 
Furthermore, we compute the dimer susceptibilities for 
different pattern candidates. While local pinwheel corre- 
lations tend to get enhanced by the RG flow, long-range 
VBC orders are rejected as seen by a strongly inhomoge- 
neous pattern response (Fig. [3]). This suggests that the 
initial RG bias for VBC flows away towards a random va- 
lence bond (RVB) liquid type state, which is consistent 
with spin liquid proposals for the KHM. 

Model. The Hamiltonian of the J\-Ji kagome Heisen- 
berg model (KHM) is given by 

^khm = j\ s^Sj + J2 s^Sj, (1) 

where (ij) and ((ij)) denote nearest neighbor and second 
neighbor pairs, respectively. We parametrize the cou- 
plings by J\ = J cos a and J 2 = J sin a which enables 
us to characterize each point in parameter space by a 
single angle a (Fig. [I]), with < a < 2ir. For a = 90°, 



the system decouples into three independent Kagome lat- 
tices which implies that the system at this point has the 
same physical properties as at a = 0. For its classical 
counterpart [30], the J1-J2 KHM exhibits four types of 
magnetic order: (i) the planar q = 0-Neel state with a 3 
site unit cell which appears for purely antiferromagnetic 
interactions at 0° < a < 90°, (ii) the non-planar cuboc 
state with a 12 site unit cell located at 90° < a < 161.6°, 
(iii) a ferromagnetic phase at 161.6° < a < 270° and (iv) 
the planar \/3 x V^-Neel state with a 9 site unit cell in 
the region 270° < a < 360° (the phases are also found in 
the quantum phase diagram Fig. [I]). These types of or- 
der correspond to different ordering-peak positions in the 
second Brillouin zone [31 j. For real-space illustrations of 
the different orders, we refer to the Refs. l30l and l32l 

Pseudofermion functional renormalization group. The 
PFFRG approach [25-28 , which we employ to obtain 
the quantum phase diagram displayed in Fig. [I] starts 
by reformulating the spin Hamiltonian in terms of a 
pseudo fermion representation of the spin- 1/2 operators 

S» = WLapfKefa =t4, » = W) with 

fermionic operators and f± and Pauli-matrices a^. 
This allows to apply Wick's theorem leading to standard 
Feynman many-body techniques. We further introduce 
an infrared frequency cutoff A in the fermionic propaga- 
tor. The FRG ansatz then formulates equations for the 
evolution of all m-particle vertex functions [33 under the 
flow of A. To reduce the infinite hierarchy of coupled 
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FIG. 2. Non-magnetic phase for \^- \ <C 1: (a) A-flow of the largest fc-component of the magnetic susceptibility for a — 0. The 
flow behavior is smooth to A = 0. (Small oscillations below A « 0.3 are due the frequency discretization.) (b), (c) fc-space 
resolved susceptibility at A = 0. (% and A are given in units of Ji.) Small maxima can be seen at the K points of the second 
Brillouin zone for J2 = in (b). These maxima vanish for at J2 = 0.017Ji in (c), with strong resemblance to [21] . 



equations to a closed set, the PFFRG only includes 2- 
particle reducible two-loop contributions [34], which still 
assures a sufficient back-feeding of the self-energy cor- 
rections to the two-particle vertex evolution. A crucial 
advantage of the PFFRG as compared to Abrikosov-type 
spin-RPA methods is that the summations include ver- 
tex corrections between all interaction channels, i.e., the 
two-particle vertex includes graphs that favor magnetic 
order and those that favor disorder in such a way that it 
treats both tendencies on an equal footing. A numerical 
solution of the PFFRG equations requires i) to discretize 
the frequency dependencies and ii) to limit the spatial 
dependence to a finite cluster. In our calculations, the 
latter typically includes a correlation area (cluster size) 
of 317 lattice sites of the kagome lattice. The onset of 
spontaneous long-range order is signaled by a sudden 
breakdown of the smooth RG flow, while the existence 
of a stable solution indicates the absence of long-range 
order [25j [26] . From the effective low-energy two particle 
vertex, we obtain the spin susceptibility with a high mo- 
mentum resolution (Fig. [I] and Fig. [2| and, in case of a 
magnetically disordered regime, track possible VBC or- 
ders [27] by computing the response functions to different 
dimer patterns (Fig. [3]). 

Quantum phase diagram. As shown in Fig. IT] we find 
ferromagnetic, cuboc, q = 0, and y/3 x Vo-Neel or- 
der in the quantum J1-J2 KHM. Fig. [I] (bottom cen- 
ter) shows the RG flow in the four magnetically or- 
dered phases. Clear instability breakdowns as resolved 
by strong kinks are seen in the ferromagnetic and in the 
y/3 x ^/3-Neel phase. In the q = 0-Neel phase and in the 
cuboc phase, however, the instability features are less 
pronounced, which is consistent with a small magneti- 
zation. We conclude that disorder fluctuations are still 
important in these regions. Fig. [I] shows susceptibility 
profiles of the magnetically ordered phases (the plots are 
taken for the cutoff value A c at the instability break- 
down). It exhibits distinct peak structures as expected 
for the different orders [31]. Note that for the ferromag- 



netic and the y/S x ^/3-Neel phase, the PFFRG also ac- 
curately resolves the expected subdominant peaks. 

Due to quantum fluctuations, magnetically disordered 
phases complement the ordered phases in the J\- J 2 KHM 
quantum phase diagram. We find a non-magnetic region 
for 338° < a < 35° (Fig. [IJ. A similar non-magnetic 
phase is found for 56° < a < 124°. In addition, we find 
indication for a small disordered phase at a « 161. En- 
hanced uncertainties near the phase boundary between 
the cuboc and the ferromagnetic phase, however, make 
it hard to resolve the range of this phase. 

The comparably small regime of unambiguous cuboc 
order is an important information in light of recent ex- 
perimental findings for Kapellasite [23 . Our results con- 
firm that the only way to accomplish such a phase in the 
quantum J\- J 2 KHM is the existence of a ferromagnetic 
J\ and a considerable antiferromagnetic J2, while even 
longer-range Heisenberg couplings might additionally be 
important. In particular, we find the cuboc phase in 
a regime with enhanced quantum fluctuations, suggest- 
ing that comparably small changes of system parameters 
should induce a significant change in Tn and the general 
magnetic susceptibility profile. The PFFRG is ideally 
suited to track the evolution of the susceptibility peaks 
in the cuboc domain, which we defer to a later point. 

Disordered phase for |^| < 1. The J x KHM (a = 0) 
has been frequently studied in the literature, and is sup- 
posed to be close to the parameter regime where Her- 
bertsmithite is located. There, we find that the RG flow 
remains stable in the entire flow regime, such that the 
RG equations can be smoothly integrated down to A = 0. 
Fig. [2|i depicts the largest susceptibility component (i.e., 
the one at the K points of the second Brillouin zone). No 
sign of any instability breakdown is seen during the flow, 
which is a signature for a non-magnetic ground state. 
This property along with the indication for a spin gap 
is also reflected in the spin-spin correlations (inset of 
Fig. [2^i). Our data show an almost perfect exponen- 
tial decay of the correlations in real space, from which 
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FIG. 3. Dimer configuration for (a) the HVBC [35] and (b) the modified HVBC before (left) and after (right) the RG flow. 
Red lines correspond to strengthened and green lines to weakened bonds, where the magnitude of the response is encoded in 
the thickness of the lines. The blue lines marks the unit cell of the pattern, (a) The local pinwheel structure gets enhanced 
during the flow, while the outer bond responses are inhomogeneous indicating the rejection of long-range order, (b) A modified 
HVBC with 60° rotation symmetry shows similar features as (a). 



we determine a small correlation length of only £ = 0.98 
lattice constants in correspondence with DMRG calcu- 
lations [12]. The full fc-space resolved susceptibility in 
Fig. [2}d reveals further information about the magnetic 
properties at a = 0. The susceptibility is mainly con- 
centrated at the boundary of the second Brillouin zone 
and is almost constant along the whole edge. Fluctua- 
tions at such large k are consistent with the small cor- 
relation length. The susceptibility along the edge shows 
only small maxima at the K points, which have also been 
found in DMRG calculations [12]. As a result, subleading 
y/3 x y/3 fluctuations are slightly preferred as compared to 
other fluctuations at the Brillouin zone boundary. Most 
of the discovered features such as a broad susceptibil- 
ity profile are shared by neutron scattering susceptibility 
measurements in Herbertsmithites, and have been inter- 
preted as signs of spin fractionalization which can emerge 
in various spin liquid scenarios [2T] , One major discrep- 
ancy, however, is the fact that the measurements do not 
find any maxima at the K points. By increasing J2, we 
discover that the previous maxima at a = vanish and 
look similar as the experimental finding: Fig. [2|c) shows 
the spin susceptibility at J2/J1 = 0.017, indicating that 
the absence of the boundary peaks in Herbertsmithite 
can be explained by the presence of a small, but finite 
J2 > coupling. 

The natural competitor for a spin liquid in such a mag- 
netically disordered phase is VBC order. For a = 0, we 
have considered generalized susceptibilities which mea- 
sure the propensity of the system to form a specific VBC. 
As has been shown by dimer expansions [35] , the hon- 
eycomb VBC (HVBC) (Fig. |3^i) is the most promising 
dimerization pattern for the J\ KHM. It is a longstanding 
question whether or not the ground state of this system is 
a HVBC. In order to calculate valence bond susceptibili- 
ties for such dimer coverings within our RG approach, we 
add a small perturbation to the model in which 
increases ( J\ —> J\ + 8) or decreases ( J\ —> J\ — S) the 
couplings on the nearest neighbor bonds, according to the 
dimer pattern. If the system supports the dimer pattern, 
the strong bonds become stronger and the weak bonds 



become weaker during the flow [31]. 

The RG flow of the HVBC pattern is depicted in 
Fig. [3^i. The left picture denotes the modified bond 
strengths at the start of the RG flow. The right picture 
shows the bond strengths after the flow at A = 0. As is 
clearly visible, the pinwheel bonds around the rotation 
center of the HVBC structure are strongly and homo- 
geneously enhanced. This has been noticed before, as 
the pinwheel is an approximate local eigenstate for the 
KHM [35]. For the complete HVBC, however, we find 
decreasing and less homogeneous response profiles as we 
consider the bonds near the boundary of the so-defined 
HVBC unit cell. As a further attempt, we have inves- 
tigated a modified HVBC pattern where the pinwheel 
structure stays unchanged, but the arrangement of the 
outer bonds is modified (Fig. [3]o) . The results we ob- 
tain show the same degree of inhomogeneous response as 
for the HVBC. Keeping the pinwheel structure, we have 
also investigated further possible VBC patterns [31] . The 
unified picture emerging from our studies is that while 
the pinwheel structures are reasonable guesses for the lo- 
cal correlation profile in the KHM, we do not find any 
sign of long-range VBC order, as we only find inhomoge- 
neous dimer responses for any VBC pattern considered. 
More so, the finding of many competing VBC pattern 
candidates frustrating each other closely resembles the 
short-range RVB liquid hypothesis, which can seed the 
emergence of spin liquid physics [36] . 

To summarize, we have obtained the quantum phase 
diagram of the J1-J2 KHM and investigated its magnetic 
and non-magnetic phases. In light of recent experiments 
on Kapellasite, it will be interesting to further investi- 
gate the interplay of order and disorder in the cuboc 
order regime. Furthermore, we have provided model ev- 
idence for explaining the features of the susceptibility 
measurements in Herbertsmithite, and specifically link 
the absence of boundary susceptibility enhancements to 
the existence of small but finite J2 coupling. We have also 
found indications for RVB liquid physics in the J\ KHM, 
rendering the related Herbertsmithite a prime candidate 
for observing a spin liquid phase in nature. 
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Supplementary material 

CLASSICAL PHASE DIAGRAM 

In the following, we discuss the classical J1-J2 KHM 
as obtained in the large-spin limit in more detail. Fig. [4] 
shows the corresponding phase diagram comprising q = 
order, cuboc order, ferromagnetic order, and y/3 x y/3 or- 
der. All phase boundaries coincide either with the J\ or 
the J 2 axis, except for the transition between the cuboc 
and the ferromagnetic phase, which lies at J 2 = — Ji/3. 
The peak structures of the Fourier transforms of these 
types of order are also shown in Fig. [4] We emphasize 
that Fourier transformed quantities on the Kagome lat- 
tice have the periodicity of the second (extended) Bril- 
louin zone, indicated by black hexagons in Fig. [4j Un- 
like the cuboc and q = order, the ferromagnetic and 
y/3 x y/S order cannot be decomposed into harmonics 
with equal wave vectors | k | . These types of order require 
two inequivalent harmonics, where the dominant (sub- 
dominant) contribution resides at the boundary of the 
second Brillouin zone (inside the second Brillouin zone). 

DIMER CORRELATIONS 

In order to investigate the properties of the magneti- 
cally disordered regimes of the quantum J1-J2 KHM we 
have probed the system with respect to its propensity to 
form various dimerized VBCs. By definition, a VBC is 
a periodic arrangement of dimer bonds on pairs of near- 
est neighbor sites (i, j), such that each site belongs to 
exactly one dimer. In the following, for a given configu- 
ration, S denotes the set of all such pairs of sites while W 
is the set of all other nearest neighbor bonds. Within our 
FRG framework, a conceptually simple way to calculate 
generalized dimer susceptibilities ^ dimer 5 measuring the 
tendency of the system to support a specific VBC, is to 
add a small perturbation to the Hamiltonian, 

(ij)es (ij)ew 

which strengthens the couplings Jij on all dimer bonds 
in S [Jij —> Jij + 5 if G S] and weakens all other 
nearest neighbor couplings [J^ — » — 6 if G W]. 
Such a modification affects the initial conditions of the 
RG flow at large cutoff scales A. As A is lowered, we keep 
track of the evolution of all static nearest neighbor spin- 
spin correlations Xij ~ 

((SiSj))^ = 0). We then define 
the dimer susceptibility for a given pair of adjacent sites 
(hj) by xtj mer = Xij-Xm where Xm is a properly chosen 
mean value over all nearest neighbor bonds (see below). 
If the absolute value |Xij mer | is small, the system tends 
to equalize (i.e. reject) the perturbation on that link, 




FIG. 4. The phase diagram of the classical J1-J2 KHM with 
q = order, cuboc order, ferromagnetic order, and y/3 x 
y/3 order. The left and right sides of the figure show the k- 
space positions of the classical ordering peaks corresponding 
to these types of order (red dots) . The black (gray) hexagons 
represent the boundaries of the second (first) Brillouin zone. 
Numbers indicate the relative peak heights of the different 
ordering peaks, which always add up to one within the second 
Brillouin zone. 

while a large value indicates that the system supports 
the dimerization. In Figs. [3] and [5J |Xij mer | is encoded 
in the thickness of the lines connecting nearest neighbor 
bonds. Furthermore, the color (green, red) indicates the 
the sign of Xij mer which distinguishes between weakened 
and strengthened bonds. 

For the HVBC introduced in the paper, as well as for 
most of the other patterns that we have tested, the unit 
cell consists of 6 inequivalent strengthened nearest neigh- 
bor bonds with G S and 18 inequivalent weakened 
nearest neighbor bonds with G W. It is convenient 
to define the mean value of the spin-spin correlations by 

Xm = \ J2 3xi i + Xi i / 36 ' ( 3 ) 

\(i,j)es (i,j)ew J 

where the sums only run over inequivalent links in a 
unit cell. In Eq. ([3| the strengthened bonds contribute 
with an additional factor of 3. This ensures that in total 
strengthened and weakened bonds are equally weighted 
such that at the beginning of the RG flow (i.e., at large 
A) |Xij mer | is constant on all bonds (left-hand pictures 
in Figs. [3] and [5|. In general, not only the A-evolution 
of the magnitude of the dimer susceptibilities |Xij mer | is 
relevant but also their spatial homogeneity: The forma- 
tion of a VBC is indicated by a large and homogeneous 
dimer response. On the other hand, an inhomogeneous 
response shows that certain bonds might have a bias to- 
wards dimerization while the dimer pattern as a whole is 
rejected and long-range order does not develop. 

In the main text, we have discussed the HVBC as well 
as a slightly modified pattern (Fig. [3]). While the central 
pinwheel structure gains considerable weight during the 
RG flow, the overall response is rather inhomogeneous. 
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FIG. 5. Further valence-bond configurations that we have 
studied. Again, red lines correspond to strengthened and 
green lines to weakened bonds. Blue lines mark the unit cell 
boundaries (note the smaller unit cell in (c)). The pattern in 
(a) is similar to the valence bond configurations in Fig. [3] i.e., 
it exhibits "perfect hexagons" at the unit cell corners, (b) 
represents a pattern of clockwise and counterclockwise rotat- 
ing pinwheels. A pattern of pinwheels with the same sense of 
rotation is shown in (c). 



Keeping the 36 site unit cell and the pinwheel, there ex- 
ists another related pattern, which we have also tested, 
see Fig. [5^i. The susceptibility distribution shows the 
same degree of inhomogeneity as in Fig. [3j Note that the 
corners of the unit cells in Figs. [3] and [5^i exhibit a local 
dimer arrangement which has been referred to as "perfect 
hexagons" (i.e. hexagons with three dimer bonds along 
its edges). Obviously these features represent the main 
source of small and inhomogeneous responses in Figs. [3] 
and [5^. 

Given these results it is interesting to study valence 
bond configurations which only consist of pinwheels and 
no "perfect hexagons". Fig. [5}d shows such a pattern 
where each strengthened dimer bond belongs to a pin- 
wheel (this is best seen by putting together various unit 
cells). Together with Figs. [3] and [5^l, these four patterns 
are the only possible valence bond configurations with 
a 36 site unit cell, a pinwheel structure and (at least) 
120° rotation symmetry. Note that as in Fig. [5)3, pin- 
wheels have a handedness (clockwise and counterclock- 
wise rotating pinwheels occur in the ratio 1:3). Again, 
while the pinwheels become the salient features during 
the RG flow, other bonds exhibit small responses, ruling 
out long-range order of this type. 

Finally, we have tested a dimer pattern which only con- 
sists of pinwheels with the same handedness, (Fig. [5^). 
This gives a unit cell of only 12 sites. The response is 
very similar to Fig. [5}d, which in total confirms the find- 
ing of a low energy scenario involving many competing 
valence bond configurations. 
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